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VO ! Abstract 

o 

' Although there is now a good measure of agreement between 

Monte Carlo and high-temperature series expansion estimates for Ising 
I (n = 1) models, published results for the critical temperature from se- 

ries expansions up to 12th order for the three-dimensional classical 
Oh. Heisenberg (n = 3) and XY (n = 2) model do not agree very well 

(— I I with recent high-precision Monte Carlo estimates. In order to clarify 

this discrepancy we have analyzed extended high-temperature series 
. expansions of the susceptibility, the second correlation moment, and 

I the second field derivative of the susceptibility, which have been de- 

' rived a few years ago by Liischer and Weisz for general 0(n) vector 

spin models on L)-dimensional hypercubic lattices up to 14t/i order in 
K = J/ksT. By analyzing these series expansions in three dimensions 
with two different methods that allow for confluent correction terms, 
we obtain good agreement with the standard field theory exponent 
estimates and with the critical temperature estimates from the new 
high-precision Monte Carlo simulations. Furthermore, for the Heisen- 
berg model we also reanalyze existing series for the susceptibility on 
the BCC lattice up to 11th order and on the FCC lattice up to 12th 
order using the same methods. 



1 Introduction 



In the past few years considerable progress has been made in developing 
very efficient Monte Carlo (MC) simulation techniques [Q. This allows high- 
precision computations of the critical coupling and the critical exponents 
of continuous phase transitions with an accuracy that is comparable with 
the widely accepted estimates derived from field theory [0, |^. The third 
and oldest approach to extract information about the critical properties of 
those systems are analyses of high-temperature series expansions. For some 
standard models (with notable exceptions including the three-dimensional 
Ising model Q and certain two-dimensional systems), however, the critical 
coupling and the critical exponents calculated by this method have much 
larger error bars and are more vulnerable to systematic errors. In order to 
improve this situation two points are important. First, more refined methods 
of analysis than in the pioneering works must be employed, and second it 
is obvious that longer series are needed. The first point should cause no 
problem anymore for continuous phase transitions since over the years many 
greatly refined methods have been developed that take into account various 
confluent correction-to-scaling terms and are now available on a routine basis 
[^, [7[. Confluent corrections-to-scaling arise from irrelevant operators and 
their neglect can bias critical coupling and critical exponent estimates. The 
generation of longer series, however, is still a very demanding numerical and 
computational problem, even though it appears to be trivial in principle. 

Significant progress in series generation has been made with star graph 
and no- free-end (NFE) graph |^ |T^, rT| enumerations which lead to medium 
length series in general dimensions for many systems. However, these ap- 
proaches are limited by the order of the existing graph table and not all 
problems have star or NFE formulations; even when these exist, the imple- 
mentation can be quite complex. For the classical 0{n) vector spin models 
an important step forward has been made by Liischer and Weisz |12|, who 
applied linked cluster expansion techniques to compute the expansion coef- 
ficients of the susceptibility, the second correlation moment and the second 
field derivative of the susceptibility on D-dimensional hypercubic lattices up 
to the lAth order in the expansion parameter K = J/ksT and provided ex- 
plicit tables for l<n<4, 2<D<4. Moreover, Butera et al. |13| observed 
that the symmetry of these models implies (Schwinger-Dyson) identities be- 
tween correlation functions that allow a recursive computation of the series 
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expansion coefficients and reveal their structure as function of n. Combining 
their result with those of Liischer and Weisz they were able to give the ex- 
pansion coefficients in general form as ratios of polynomials in n. Although 
still one term shorter than the NFE tables |1^ , and three terms below the 
star graph series of Singh and Chakravarty P], these methods can be used to 
generate longer series directly, requiring only larger computer memory and 
not preexisting graph tables. 

The motivation to analyze the extended high-temperature series expan- 
sions of the Heisenberg model comes from two recent high-precision MC 
simulation studies [|1^, [1^ of this model on simple cubic (SC) lattices which 
gave significantly larger values for the critical coupling than previous esti- 
mates based on analyses of series expansions up to 12th order [16-19], and 
transfer matrix MC studies see Table 1 (also included is newer MC data 
[ PH] , that was obtained after completion of our work). There are two sources 
for the expected improvement. First, on hypercubic lattices two more terms 
of the series are known and second, more refined methods taking into account 
confluent correction terms are available. For the latter reason we also rean- 
alyze the long-known but shorter series for the susceptibility on the body 
centered cubic (BCC) and face centered cubic (FCC) lattices. Finally, we 
present analyses of the new longer series for the XY {n = 2) model on the 
SC lattice. 



2 Model and observables 



We consider the classical 0{n) symmetric Heisenberg model with partition 
function 



n 



exp 



Si ' S j 



(id) 



(1) 



where K = J/ksT is the reduced inverse temperature, denotes nearest- 
neighbor pairs, and Qi is the surface of the n-dimensional unit sphere as- 
sociated with the degrees of freedom of the n-dimensional unit spins Si at 
each site of a regular three-dimensional lattice. In this paper we investigate 
the new longer series for the Heisenberg (n = 3) model on a SC lattice, and 
reanalyze existing series for the BCC and FCC lattices. Further we also 
study the new longer series for the XY {n = 2) model on a SC lattice. In or- 
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der to estimate the critical couplings and exponents we concentrate on three 
observables, the susceptibility 

2 



X 



) 



1 + a^t^^ +b^t + ... 



(2) 



the second correlation moment 

i 



X- 



, T,ii'^{so ■ Sj) 



1 + a-m{2)t 



(3) 



and the second field derivative of the susceptibility 
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X 



(4) 



n(n + 2' 



^(So ■ Si Sj ■ Sk)c 



where (...) denotes expectation values with respect to the partition function 
(|l]) and the subscript c in (|) stands for the connected part. The second lines 
in (D) - (§) give the assumed critical behavior where t = — > is the 
distance from the critical point in the high-temperature phase, 7, u and f3 are 
the standard critical exponents of the susceptibility, correlation length and 
magnetization, respectively, and the terms in square brackets describe the 
leading confiuent and analytic correction terms. In (Q) we have made use of 
the relation A = 7+/?, where A is the gap exponent. In the high-temperature 
phase these observables can be expanded as 



m 



(2), 



r=l 

00 



n,K) = J2br{n)K' 



r=l 



X^'\n,K) 



-2 + Y,dr{n)K' 



r=l 



(5) 
(6) 
(7) 



n{n + 2) 

defining the coefficients ar{n), hr{n) and dr{n) computed in refs. [0, |13 
For the convenience of the reader we have compiled their numerical values 
for n = 2 and = 3 in Tables 2 and 3. 
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3 Methods of analysis 



We analyze the series given in Tables 2 and 3 with two different methods ^2 
that allow for confluent and analytic correction terms. Taking the suscep- 
tibility as a generic example (and suppressing subscripts) we thus assume a 
critical behavior of the form 

X = At-^ [l + at^' +bt + (8) 

where Ai = uuj 0.55) is the confluent correction exponent and bt a (sub- 
leading) analytic correction term. The non-universal amplitudes A, a, b are 
assumed to be constant. The . . . inside the brackets indicate further higher 
order corrections of the form t^'", which we neglect in our analysis. 

In the method referred to as Ml, first the leading singularity is removed 
by forming 



Ar 



Aiat^' +bt + ... . (9) 



Then Pade approximants are applied to the logarithmic derivative of B, 

d\nB _ Ai(7 - Ai)at^i-i + (7 - 1)6 
dt ~ t(Aiat^i-i + b) ' 



(10) 



yielding for given Kc the confluent correction exponent Ai as function of 7, 
Ai = Ai(7). The optimal set of values for the parameters 7 and Ai is 
determined visually from the best clustering of different Pade approximants. 

In the second method referred to as M2, Pade approximants in a new 
variable 

^ = 1 _ (1 _ K/K,)^^ = 1 - {t/K,t^ (11) 

are applied to 

(91nx Aiat^i+6t 
t = -J + 



dt 1 + at^i + bt 

^ l-ir,a(y-l)-ir,(y-l)VA.' ^ ) 
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yielding for given Kc the exponent 7 as function of Ai, 7 = 7(Ai). Again 
the clustering of different Pade approximants is used to select the optimal 
set of parameters. 



The two methods are complementary and as stressed in App. D of ref. [22 



should always be used in conjunction to avoid spurious results due to so-called 
resonances at values of Ai/n, n = 2,3, .. . in the otherwise more accurate 
method M2. The analysis was made with the help of the recently developed 
VGS program package that makes extensive use of the graphic features 
of an X-window environment and allows easy and efficient scanning of the 
three-dimensional parameter space. 



4 Results 

4.1 Heisenberg {n = 3) model 

SC lattice: As mentioned in the introduction our main emphasis was on 
the Heisenberg model on a SC lattice since recent high-precision MC sim- 
ulation studies [0, |T5[ were at odds with previous high-temperature series 
expansion analyses [16-19]. In particular the critical coupling Kc turned out 
to be significant larger than widely accepted series estimates based on ex- 
pansions up to 12th order; see Table 1. Our main result from analyses of the 
longer 14 terms series using methods Ml, M2 is that we can clearly confirm 
the MC estimates of Kc- More precisely for all three series we get consis- 
tent results from methods Ml and M2, and the three estimates for Kc vary 
only weakly: Kc = 0.6928 from analyses of -^c = 0.6930 from m*^^) and 
Kc = 0.6928 from x^^^- Taking the average of these three values as the final 
result we get 

Kc = 0.6929 ± 0.0001 (SC lattice). (13) 

To illustrate the method of analysis we show for the susceptibility in Fig. 1 
graphs of the highest near diagonal Pade approximants to the critical expo- 
nent 7 in the three-parameter space Kc, Ai, 7 computed according to method 
M2. A two-dimensional plot of the central slice at Kc = 0.6928 is shown in 
Fig. 2(b). The corresponding plot for method Ml is displayed in Fig. 2(a). 
From the point of best clustering of the different Pade approximants shown 
in Fig. 2 we read off 

7 = 1.400 ±0.010, (14) 
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and Ai = 0.7 ± 0.2. Similar analyses of the series for m*^^-* yield 7 + 2z/ 
2.825 ± 0.020 or inserting 



= 0.712 ±0.010, (15) 

and from x^^^ we get 87 + 2/3 = 4.925 ± 0.020 or using ([l|) 

/3 = 0.363 ± 0.010. (16) 

Using the scaling relation a + 2/3 + 7 = 2 and the estimates (p!^,(|T6|) we 
calculate 

a = -0.125 ±0.020. (17) 

Since we have three independent estimates of critical exponents this result 
can be used to test the hyperscaling relation a = 2 — Du. Using the estimate 
(|T5|) we obtain 

a = -0.136 ±0.030, (18) 

in good agreement with (|l^), thus supporting the hyperscahng hypothesis. 
Similarly, the scaling relation 5 = 1 ± 7//? gives 

(5 = 4.86 ±0.10, (19) 

while the hyperscaling relation 7/1/ = 2 — t] = D{6 — l)/{6 + 1) yields a 
comparison between central estimates of 7/// = 1.966 from the values quoted 
above and the r.h.s. of the scaling relation 

D(5- l)/(5± 1) = 1.975, (20) 

again in good agreement with each other. Our results for the critical expo- 
nents are summarized in Table 4, where they are compared with the standard 
field theory values and the results from recent MC simulations. 

BCC lattice: The susceptibility series consists of only 11 terms, but 



the overall behavior is similar. We find optimal convergence at 

= 0.4867 ± 0.0001 (BCC lattice), (21) 

again with 7 ^ 1.4 but with a lower correction-to-scaling exponent Ai than 
was seen in the SC case. We quote central estimates of Ai ^ 0.6 from Ml 
and Ai ^ 0.5 from M2. 
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FCC lattice: For the FCC lattice the 12th order susceptibility series was 
analyzed, including corrections-to-scaling, in ref. |T^. It was found that the 
amplitude of the confluent correction (with Ai = 0.55 held fixed at the RG 
value [0]) was very small, and that the analytic correction was the dominant 
one. We find 

= 0.31475 ± 0.00010 (FCC lattice) (22) 



and 7 ~ 1.39, in good agreement with |T9[. This 7 is a little lower than our 



values on the other lattices, and closer to the values of other calculations. In 



contrast to |19|, we saw clear evidence of a non-analytic correction-to-scaling 



at Ai ^ 0.6 from the M2 study of a first derivative of the susceptibilty series. 



4.2 XY {n = 2) model 

SC lattice: For the XY model we have only analyzed the new longer series 
for the simple cubic lattice. In this case the series for the susceptibility 
turned out to be not well-behaved and it was very difficult to get precise 
estimates of the critical parameters. With this caveat in mind we estimate 
Kc = 0.45407 and 7 = 1.325. On the other hand the series for m^'^^ and 
X^^^ behaved similar to the Heisenberg model, i.e., both methods Ml and M2 
gave consistent results and the estimates of Kc from both series agreed with 
each other, 

Kc = 0.45414 ± 0.00007 (SC lattice). (23) 

While the previous estimate = 0.4539 from series analyses is again 
lower (and clearly below the error limits of the present study), our value (p^) 
is consistent with recent Monte Carlo studies which gave = 0.45421(8) p4| 
using multiple and Kc = 0.4542(1) using single cluster simulations. For 
the exponents we obtain central estimates of 7+2z/ = 2.67 from the expansion 
of m^^\ and 37 + 2/5 = 4.67 from the expansion of x^^^- The exponents 
calculated from these estimates are u = 0.673, 7/2^ = 1.970 = 2 — 1], and 
jS = 0.348. These values are again consistent with field theoretical estimates 
[0, The scaling relations yield a = —0.020 and 6 = 4.81. The hyperscaling 
relations result in a = —0.018, and ^jv = D{6 — l)/{6 + 1) = 1.968, again 
in good agreement with our previous values. 
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5 Concluding remarks 



Analyzing new longer series for the Heisenberg {n — 3) model using more 
refined methods than in early works we obtain for the SC lattice critical 
parameters that are in good agreement with completely independent results 
from two recent MC simulations. Our reanalyses of existing series for the 
FCC and BCC lattice indicate that the improvement comes mainly from the 
refined methods that are able to take into account confiuent correction terms. 
With 14 (or even only 12 or 11) terms these series are, however, still too 
short to compete with the accuracy achieved by field theoretical methods for 
critical exponents, or with the precision claimed from simulations. However, 
the results clearly show that there remain no major discrepancies between 
series estimates and other calculations. Longer series clearly stabihze and 
thus increase the reliability of the estimates along the lines discussed here, 
and it therefore would be very desirable to have a few more terms available. 
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Tables 





method 


authors 


0.692 


8 terms HTS 


Wood and Rushbrooke (1966) |16| 


0.692(4) 


9 terms HTS 


Joyce and Bowers (1966) [|1^] 


0.6916(2) 


9 terms HTS 


Ritchie and Fisher (1972) p| 


0.6924(2) 


12 terms HTS (Fade) 


McKenzie et al. (1982) |g 


0.6925(1) 


12 terms HTS (ratio) 




0.6922(2) 


TMMC {n > 5) 


Nightingale and Blote (1988) ^ 


0.6925(3) 


TMMC \n > 6) 




0.6929(1) 


Metropolis MC 


Feczak et al. (1991) IH 


0.6930(1) 


1 Cluster MC 


Holm and Janke (1992) ^ 


0.693035(37) 


multiple 1 Cluster MC 


Chen et al. (1993) fH 


0.6929(1) 


14 terms HTS 


this work 



Table 1: Estimates of the critical coupling of the Heisenberg (n = 3) 
model on a simple cubic lattice from various sources (HTS: high-temperature 
series analysis, TMMC: transfer-matrix Monte Carlo, MC: Monte Carlo sim- 
ulation). 
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order r 




br 




1 


3 


0000000000 


3 


0000000000 


-24 


0000000000 


2 


7 


5000000000 


18 


0000000000 


-160 


5000000000 


3 


18 


3750000000 


72 


3750000000 


-822 


0000000000 


4 


43 


5000000000 


247 


5000000000 


-3576 


8125000000 


5 


102 


3437500000 


770 


5937500000 


-13971 


7500000000 


6 


237 


0546875000 


2261 


3437500000 


-50454 


9648437500 


7 


546 


9462890625 


6360 


6650390625 


-171739 


3593750000 


8 


1252 


0048828125 


17343 


7773437500 


-557978 


9429687500 


9 


2858 


8175292969 


46158 


4210449219 


-1746304 


9972656250 


10 


6496 


1514078776 


120515 


3193033854 


-5299323 


3505303277 


11 


14735 


3746412489 


309746 


4250318739 


-15671446 


8761067708 


12 


33314 


7537746853 


785831 


2964274089 


-45336965 


5964835394 


13 


75222 


2566392081 


1971809 


9920579093 


-128702556 


1244287884 


14 


169444 


4882359232 


4901417 


5916496216 


-359396456 


8541712222 



Table 2: Expansion coefficients for tlie XY (n — 2) model high-temperature 
series for a simple cubic lattice. Given are the expansion coefficients of 

the susceptibility the expansion coefficients br of the second correlation 
moment m*^^\ and the expansion coefficients dj. of the second field derivative 
of the susceptibility x^^^ up to lAth order (for details compare text). 
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order r 




br 




1 


2 


0000000000 


2 


0000000000 


-16 


0000000000 


2 


3 


3333333333 


8 


0000000000 


-71 


7333333333 


3 


5 


4222222222 


21 


4222222222 


-246 


0444444444 


4 


8 


5185185185 


48 


7111111111 


-716 


4486772487 


5 


13 


2670194004 


100 


7336860670 


-1870 


2019047619 


6 


20 


3359905938 


196 


1285831864 


-4508 


3329617872 


7 


30 


9989637468 


365 


7050425240 


-10232 


2542817950 


8 


46 


8673402574 


660 


4991803514 


-22145 


7412271162 


9 


70 


6067866595 


1163 


5584276550 


-46128 


4203352476 


10 


105 


8320214871 


2009 


6315902889 


-93088 


6148720584 


11 


158 


2324753396 


3414 


9732182123 


-182932 


5061463846 


12 


235 


7598652836 


5725 


3717946474 


-351440 


3272602895 


13 


350 


6189575427 


9489 


5939248535 


-662121 


9818887996 


14 


520 


1310140421 


15575 


4527177723 


-1226410 


1925173962 



Table 3: Expansion coefficients for the classical Heisenberg {n ~ 3) model 
high-temperature series for the simple cubic lattice. Given are the expansion 
coefficients of the susceptibility x, the expansion coefficients br of the 
second correlation moment m^'^\ and the expansion coefficients dr of the 
second field derivative of the susceptibihty x^^^ up to lAth order (for details 
compare text). 
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method 


V 


7 


/3 


a 


b 


5f-expansion p| 


0.705(3) 


1.386(4) 


0.3645(25) 


-0.115(9) 


4.802(37) 


e-expansion 


0.710(7) 


1.390(10) 


0.368(4) 


-0.130(21) 


4.777(70) 


MC H 


0.706(9) 


1.390(23) 


0.364(7) 


-0.118(27) 


4.819(36) 


Mc m 


0.704(6) 


1.388(14) 


0.362(4) 


-0.112(18) 


4.837(11) 


MC [21] 


0.7036(23) 


1.3896(70) 


0.3616(31) 


-0.1108(69) 




this work 


0.712(10) 


1.400(10) 


0.363(10) 


-0.136(30) 


4.86(10) 



Table 4: Critical exponents for the three-dimensional classical Heisenberg 
(n = 3) model from various sources. 
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Figure Headings 



Fig. 1: Graphs of highest near diagonal Pade approximants to 7 in the 
three-parameter space K^., Ai, 7 for method M2. A two-dimensional 
plot of the central slice at = 0.6928 is shown in Fig. 2(b). 

Fig. 2: Graphs of highest near diagonal Pade approximants to 7 plotted 
against Ai at fixed K^. — 0.6928 for (a) method Ml and (b) method 
M2. 
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